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Abstract 

The superfluid density is calculated theoretically for incompressible vortex lattices in two dimen- 
sions that have isolated dislocations quenched in by a random arrangement of pinned vortices. The 
latter are assumed to be sparse and to be fixed to material defects. It is shown that the pinned 
vortices act to confine a single dislocation of the vortex lattice along its glide plane. Plastic creep 
of the two-dimensional vortex lattice is thereby impeded, and macroscopic phase coherence results 
at low temperature in the limit of a dilute concentration of quenched-in dislocations. 
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I. INTRODUCTION 



Consider either a bismuth-based, a mercury-based, or a thallium-based high-temperature 
superconductor in high enough external magnetic field so that magnetic flux lines appear, 
and so that these overlap considerably. Such materials are extreme type-II layered super- 
conductors, with mass-anisotropy ratios between the layer direction and the perpendicular 
direction of order 10~ 3 or less 1 . To a first approximation then, it becomes valid to neglect 
the coupling of magnetic screening currents between layers, as well as the Josephson effect 
between them. The description of the initial physical situation is thereby reduced to a stack 
of isolated vortex lattices or vortex liquids within each layer. 

High-temperature superconductors like those just mentioned typically also have crys- 
talline defects and inhomogeneities that act as pinning centers for vortex lines in the mixed 
phase. Within the preceding approximation, theoretical and numerical studies indicate that 
any net concentration of randomly pinned vortices results in a net concentration of unbound 
dislocations quenched into the two-dimensional (2D) vortex lattices of each layep^ 1 ^. De- 
fective vortex matter is then left as the only possible solid state of the mixed phase in 
two dimensions. It is useful to divide the last group into two classes: (i) configurations 
of vortices that contain no unbound disclinations, and (ii) configurations of vortices that 
contain some concentration of unbound disclinations 6 . The amorphous vortex glass charac- 
terized by macroscopic phase coherence falls into the second class 7 . It is believed to exist 
only at zero temperature in two dimensions, however. Defective vortex lattices with no un- 
bound disclinations, but with isolated dislocations^, or with dislocations arranged into grain 
boundaries^ ' 10 > 11 , are then perhaps left as the only solid states of the mixed phase that are 
possible in two dimensions above zero temperature. 

In this paper, we demonstrate theoretically that defective vortex lattices in two dimen- 
sions show macroscopic phase coherence in the extreme type-II limit, in the regime of weak 
random pinning. We find, in particular, that the 2D vortex lattice exhibits a net superfluid 
density if it is void of disclinations, and if only a small number of isolated dislocations are 
quenched in in comparison to the total number of pinned vortices. This result is achieved in 
three steps. First, we demonstrate in section II that a network of pinned vortices confines 
the motion of a single dislocation along its glide plane. This guarantees that the 2D vortex 
lattice remains elastic in the limit of a dilute concentration of such unbound dislocations. 
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Next, the uniformly frustrated XY model for the 2D vortex lattice is introduced in section 
III, through which a useful expression for the superfluid density is derived in terms of glide 
by unbound dislocations^. Interactions among the dislocations are notably ignored here. 
These results are then assembled in section IV, where the final formula [Eq. (|27|) ] for the 
superfluid density of the defective vortex lattice is obtained as a function of the ratio of the 
number of unbound dislocations to the number of pinned vortices. 

II. COLLECTIVE PINNING OF ONE DISLOCATION 

It is strongly believed that the vortex lattice in two-dimensions is unstable to the prolifer- 
ation of dislocations in the presence of an arbitrarily weak field of random pinning centers^. 
Let us assume this to be the case. Let us also assume that the interaction between dislo- 
cations can be neglected due to the screening action by the random pins 12 . This requires a 
dilute concentration of dislocations in comparison to the concentration of pinned vortices. 
Consider then a single dislocation in the 2D vortex lattice at zero temperature in the extreme 
type-II limit. The latter implies that the vortex lattice is incompressible. The dislocation 
can therefore slide along its glide planed, but it cannot climb across it. This would require 
the creation or the destruction of vortices, which is prohibited by the extreme type-II limit. 
Below, we shall demonstrate how randomly pinned vortices in the 2D vortex lattice act to 
pin the dislocation itself along its glide plane. 

Consider a single dislocation that can move along its glide plane in the 2D vortex lattice 
with randomly located material defects present. Assume that a small fraction of the vortices 
are localized at some subset of the pinning centers. The former is guaranteed at zero- 
temperature for a sparse array of random pinning centers compared to the density of vortices. 
A vortex that lies at a point R in the case of the perfect triangular vortex lattice will in 
general be displaced to a position R + u(R) by the action of thermal fluctuations and of the 
random pinning centers. We shall now make the approximation that the pinned vortices are 
fixed: 

u(Ri)=Vi for i = 1,2, JV pin , (1) 

where Ri is the home site of the vortex pinned down at Ri + Vi, and where N pin denotes 
the total number of pinned vortices. This approximation is valid for physics at large length 
scales compared to the effective radius of a pinning center. It then requires low magnetic 
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fields compared to the upper critical one if the radius of a pinning center is of order the 
coherence length. The energy of the pinned vortex lattice is then given by 

E=-^J d 2 R (V xu) 2 + J d 2 R X-[u-v\ (2) 

in the continuum limit, where /xq denotes the shear modulus of the unpinned vortex 
latticed, and where X(R) = E?=i \6^(R - Ri) is the field of Lagrange multipliers that 
is introduced in order to enforce each of the Ap; n constraints (0). Also, the vortex lattice 
is incompressible in the extreme type-II limit, and this requires that the displacement field 
satisfy the constraint 

V-m = (3) 

everywhere. By Eq. (J21), the equilibrium configuration uq of the dislocation confined to its 
glide plane then satisfies the field equation 

/i V x V x u + A = (4) 

everywhere. It can be used to show that the elastic energy (j2J) for a fluctuation about 
equilibrium, u = u + 5u, takes the form 

E = -fi J d 2 R (V x M ) 2 + i/i J d 2 R (V x 5u) 2 + J d 2 R SX ■ 5u, (5) 

where 5X(R) = Y^i SXi 5^ (R — Ri) is the field of the fluctuation in the Lagrange multi- 
pliers, 5Xi = Aj — Xj . 

To proceed further, it is convenient to decompose the displacement of vortices into pure 
wave and pure defect components: u = u wv + Udf- Suppose now that the dislocation is 
displaced by SRdf along its glide plane with respect to its equilibrium position. Notice then 
that the fluctuation in the defect component corresponds to a pair of dislocations with equal 
and opposite Burgers vectors oriented along the glide plane (see fig. HJ): 

8u di (R) = u d °?(R - 5R di ) - 4?^), (6) 

where uf^R) denotes the displacement field of the pure dislocation at its home site. At this 
stage it becomes important to observe that the pure wave and the pure defect components 
do not interact elastically^: / d 2 R (V x u) 2 = J d 2 R(V x u wv ) 2 + J d 2 R(V x «df) 2 . Application 
of this fact to Eq. (j3J) then ultimately yields the form 

E = E di + + i/xo / d 2 R (V x 5u wv ) 2 + J d 2 R5X ■ (5u wv + 8u di ), (7) 



for the elastic energy, where E^J = (/io/2) / d 2 R{V x u^l) 2 is the wave contribution to the 
elastic energy at equilibrium. Also, E d f = (/io/2) / d 2 R(V x «df) 2 is the elastic energy of 
the displaced pure dislocation, which is constant. 

The energy (JJJ) of the displaced dislocation is therefore optimized by minimization with 
respect to the pure wave component Su wv along with the constraints 

5u wv (Ri) = -5u d{ (Ri) for i = 1,2, N pin . (8) 

Its solution can be obtained by a straight forward generalization of the solution for a pinned 
elastic string (see Appendix). This yields 

iVpin TVpin 

5u wv (R) = - E E G±(R - Ri) ■ G-] ■ Su di (Rj), (9) 

i=l j=l 

where 

G±(R) =Y f ^Xz x q){qL)~ 2 {z X q) (10) 

is the transverse Greens function over an L x L square region with periodic boundary 
conditions, and where G~j is the inverse of the N pin x iVp in matrix G±(Ri —Rj). Notice that 
(JHJ) manifestly satisfies the constraints (jHJ) and the incompressibility requirement V-5u wv = 0. 
Also, direct substitution of the solution © into Eq. (JJJ) yields a change in the elastic energy 
due to the displacement of the dislocation equal to 

-. N pin iV pin 

5E pin = -/ioEE 6u6f(Ri) ■ Gi] ■ Su^Rj). (11) 

Z i=l j=l 

Yet Gl~j is the inverse of the 2D Greens function, G = — V~ 2 , projected onto transverse 
displacements (J5J) and onto the sites of the pinned vortices, {Ri}. If these sites are exten- 
sive and homogeneous, then they resolve unity at long wavelength: Yli=i — 1- We 
therefore have that G~J = (z|P£ 1 (— V 2 )Pi|j) at longwavelength, where P± denotes the pro- 
jection operator for transverse displacements (jHJ). Substitution into Eq. (fTTj) then yields the 
expression 

SE pin = J' d 2 R(y x 5u di ) 2 (12) 

for the change in the elastic energy due to the displacement of the dislocation, where the 
prime notation signals that the integral has an ultraviolet cut-off (R P i n ) of order the average 
spacing between pinned vortices [see Eq. (jA.7|) in the Appendix and ref. 15 ]. We conclude 
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that the displacement of the dislocation along its glide plane generates shear stress on the 
vortex lattice via the array of pinned vortices (see fig. Q). This then results in a restoring 
Peach-Kohler force on the displaced dislocation^ 3 -. 

Let us now compute the effective spring constant of the Peach-Kohler force experienced by 
the dislocation at small displacements from equilibrium due to the array of pinned vortices: 

5E pin = h pin (5R di ) 2 for n pin (5R di ) 2 < 1, (13) 

where n P i n denotes the density of pinned vortices per layer. The relative displacement field 
(jHJ) then corresponds to that of a pure dislocation pair of extent SR d { that is oriented along 
its glide plane. Without loss of generality, it is sufficient to consider a pair of dislocations 
centered at the origin, with the glide plane located along the x-axis. The displacement field 
is then given asymptotically b y 13 ' 16 

5u di (R) S (6/7T)(5 J R df )(XF/ J R 4 )/2, (14) 

where b is one of the equal and opposite Burgers vectors oriented parallel to the glide vector 
(5Rdf)x. This expression is valid in the limit of small displacements relative to the spacing 
between pinned vortices: n vin (6R d f) 2 <C 1. Substitution into expression (JT2j) for the change 
in the elastic energy yields the result 

fcpin = (2n pill 7iR 2 pin )~ 1 (n b 2 )n pin (15) 

for the effective spring constant of the Peach-Kohler force (|13|). where R pin denotes the 
natural ultraviolet cutoff of the array of pinned vortices^: n pin ■ TiR pin ~ 1. Equations fTBl 
and (|15|) represent the final result of this section. It indicates that the incompressible vortex 
lattice confined to two dimensions does not respond plastically to small shear stress^ 3 - when a 
dilute enough concentration of unbound dislocations are quenched in. Instead, the response 
to small shear stress should remain elastic, like in the pristine case^ 4 ^, due to the pinning 
of the quenched-in dislocations. 

III. UNIFORMLY FRUSTRATED XY MODEL 

The minimal description of the mixed phase in a layered superconductor is given by a 
stack of isolated XY models with uniform frustration over the square lattice^. Both the 
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effects of magnetic screening and of Josephson coupling between layers are neglected in this 
approximation. The thermodynamics of each layer is then determined by the superfluid 
kinetic energy 

E xy = ~ E E J»™s[A^ - A,}\ r , (16) 

fi=x,y r 

which is a functional of the superconducting phase 0(r) over the square lattice. The local 
phase rigidities within layers, J x and J y , are assumed to be constant over most of the 
nearest-neighbor links, with the exception of those links in the vicinity of a pinning site. 
The vector potential = (0,2nfx/a) represents the magnetic induction, B± = § f/a 2 , 
oriented perpendicular to each layer. Here a denotes the square lattice constant for each 
layer, which is of order the zero-temperature coherence length. Also, $o denotes the flux 
quantum, and / denotes the concentration of planar vortices per site. After taking the 
Villain approximation, which is generally valid at low temperature^, a series of standard 
manipulations then lead to a Coulomb gas ensemble with pinning centers that describes 
the vortex degrees of freedom on the dual square lattice^. The ensemble for each layer is 
weighted by the Boltzmann distribution set by the energy functional 

E vx = (2tt) 2 g 5Q J G® 5Q' + E V p - m \Q\ 2 , (17) 

(R,R') R 

written in terms of the integer vorticity field Q(R) over the sites R of the dual lattice in that 
layer, and of the fluctuation SQ = Q — f '. A logarithmic interaction, = — V -2 , exists 
between the vortices, with a strength J equal to the gaussian phase rigidity. Last, Vp in (-R) 
is the resulting pinning potential^ 9 -. 

The 2D Coulomb gas ensemble (fTTj) can be used to test for the presence or the absence 
of superconductivity. In particular, the macroscopic phase rigidity parallel to the layers is 
given by one over its dielectric constant: 

pf D) /J = 1 - Iim(27r/77 SW ) (SQ^Q^/kWAfw . (18) 

k — >0 

Here 5Q^ = — {Q%} is the fluctuation in the Fourier transform of the vorticity: = 
J2r Q{R)e %k ' R . Also, ?7 SW = /cbT/27tJ is the spin-wave component of the phase-correlation 
exponent, and M\\ denotes the number of points in the square-lattice grid. Now suppose that 
a given vortex is displaced by 5u with respect to its equilibrium location at zero temperature, 
uq. Conservation of vorticity dictates that the fluctuation in the vortex number is given by 
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SQ = — V • Su. Substitution into Eq. (fTS|) then yields the result 8 

pf^/Jo = 1 - (rHJrhw) (19) 
for the phase rigidity in terms of the vortex component of the phase-correlation exponent, 

R 

The latter monitors fluctuations of the center of mass of the vortex lattice^. Above, N vx 
denotes the number of vortices, while a vx = aj f 1 ^ 2 is equal to the square root of the area per 
vortex. Also, the prime notation above signals that the summation is restricted to the vortex 
lattice. To proceed further, we again express the displacement field as a superposition of 
pure wave and of pure defect components of the triangular vortex lattice^: Su = 5u wv + 5udf- 
Observe now that 8u wv = under periodic boundary conditions if rigid translations of 
the 2D vortex lattice are not possible. The latter is achieved by the array of pinned vortices 
(P) through the elastic forces (j2j). By Eq. (J20|) . we therefore have the result 

V' vx = *([ J2'su di f)/N vx al x (21) 

R 

for the fluctuation in the center of mass of the 2D vortex lattice. The degree of phase 
coherence in the pinned vortex lattice is therefore insensitive to its pure wave contribution. 

Consider now the hexatic vortex glas d 21 ' 22 , with a collection of randomly located un- 
bound dislocations that are quenched in by the random array of pinned vortices^. Suppose 
also that the temperature is low enough so that the thermal excitation of pairs of dislocations 
in the vortex lattice can be neglected. Within the elastic medium description (J2J), the pure 
defect component of the net displacement field is just a simple sum of the displacements due 
to each individual dislocation. And by analogy with the hexatic liquid phase of the pure 
2D vortex lattice^, we shall assume that interactions in between the unbound dislocations 
can be neglected, be they direct or be they transmitted through the field of pinned vortices. 
Expression (|2~T|) for the fluctuation of the center of mass of the 2D vortex lattice then reduces 

to 

vL = *([ EW] 2 )"*' (22) 

R 

where 5u d{ (R) denotes the fluctuation field of a given dislocation displaced along its glide 
plane (JHj), where = Nd{/N vx al x is the density of unbound dislocations per layer, and 
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where the overbar notation denotes a bulk average. A lone dislocation roams along its 
glide plane to an extent that is vanishingly small, however, in the zero-temperature limit: 
5Rdf — > as T — > 0. Without loss of generality, we can then use the asymptotic expression 
(|14|) for the corresponding fluctuation in the displacement field, 5u^ . The fluctuation in 
the center of mass of the 2D vortex lattice (|22j) is dominated by the "diagonal" on-site 
contribution, ([X^ ' |<^4f^l 2 )) which yields the estimate^ 

V ' vx = n d{ (\5R di \ 2 )(b/2a vx ) 2 \n R /a di . (23) 

Here a d { is the core diameter of a dislocation, while Rq is an infrared cut-off. The above 
logarithmic divergence associated with the latter scale justifies the neglect of the contribu- 
tion to the fluctuation in the center of mass (}2"2")) by the autocorrelator (8v$ ■ 8u^ '} at 
different points. This is due to the fact (i) that — (<5t% (a) • 5i% (&)) must decay faster than 
[(\5utf\ 2 )/27cln(R /a d{ )](a vx /R ab ) 2 by Eq. (J22|) because r]' vx > 0, and to the fact (m) that 
the former autocorrelator is short range as a result of disordering by the quenched-in dislo- 
cations. Last, given that expression (}2"3*|) was obtained by neglecting interactions in between 
isolated dislocations, it is natural to assume that the infrared scale Rq that appears there is 
set by their density, n d f- 

IV. LOW-TEMPERATURE PHASE COHERENCE 

We shall now assemble the results of the previous two sections and compute the macro- 
scopic phase rigidity of a defective vortex lattice in two dimensions. Recall that the en- 
ergy cost for a small displacement of the dislocation along its glide plane takes the form 
5E d f = hk P i a \5R d i | 2 , where k pin is the effective spring constant (JT5J) due to the randomly 
pinned vortices. This approximation for the elastic energy of the dislocation along its glide 
plane is valid in the zero-temperature limit, T — > 0, where it amounts to a saddle-point 
approximation for the Boltzmann weight in the thermal average (\5Rdf\ 2 )- The periodic 
Peierls-Nabarro potential energy along the glide plane of the dislocation shall be neglected 
for the moment^. Application of the equipartition theorem to expression (|2Hjl for the fluc- 
tuation of the center of mass of the 2D vortex lattice then yields the result 

V ' vx S (k B T)(n di /k pin )(b/2a vx ) 2 \nR /a di , (24) 



9 



which notably vanishes linearly with temperature. Substitution into expression (|19|) in turn 
yields the result 

p {2D) /J = 1 - (27rJ )(n df /A; pin )(V2a vx ) 2 ln/?o/«df (25) 

for the 2D phase rigidity in the zero-temperature limit. Substitution of the result (|15[) for 
the effective spring constant above yields the final formula for the macroscopic 2D phase 
rigidity near zero temperature, 

P [ s D) ~, f R 2 ,irJ N di (R \ 

—— = 1 - n pin • ttR in r — In — , 26 

Jo Ato< x A/pin \ a df/ 

in terms of the number of unbound dislocations, iVdf, the number of vortices, iV vx , and 
of the number of pinned vortices, iVp in , in an isolated 2D vortex lattice. Recall now the 
estimate for the shear modulus of the unpinned vortex lattice in the extreme type-II limn—^, 
/ioflvx = (n/tyJo- Substitution into expression (|26|) for the superfluid density then yields the 
yet simpler result 

Pf D) ^ ^ .N^JR^ 2 



l-(2n pin -7r^ i J-^-ln near T = 0. (27) 

^0 iV pin V <2df / 



The effect of the Peierls-Nabarro potential energy with period b that any dislocation expe- 
riences along its glide plane has been neglected above^. It becomes useful in this instance 
to define the temperature scale 

k B T = k pin b 2 ~ (fi b 2 )(n pin b 2 ), (28) 

at which point thermally induced excursions of the dislocation about its home site are typi- 
cally of the size of a Burger's vector, b. Notice first that fceTo is an extremely small fraction 
of the elastic energy scale [i§b 2 if the concentration of pinned vortices is dilute. Typical 
excursions of the dislocation will be large compared to the Burgers vector at temperatures 
above this extremely low scale: on average, |5i?df| > b at T > T . In such case, the periodic 
Peierls-Nabarro potential can be neglected because the thermal motion of the dislocation 
becomes insensitive to its relatively short period b. The periodic potential will take effect at 
extremely low temperature T < T , on the other hand, in which case it will tend to localize 
the dislocation even further about its home site. 

Let us finally close the chain of calculations by estimating the ratio of the number of 
unbound dislocations that are quenched into the vortex lattice to the number of pinned 
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vortices. Notice that it determines the phase rigidity (|27|) of the defective vortex lattice at 
T > Tq. A variational calculation by Mullock and Evetts finds that this ratio is given by 



where f pm denotes the maximum pinning force, and where a' d{ is of order the core diameter of 
a dislocation in the vortex lattice^. This result is valid only in the collective-pinning regime 
at iVdf < iVp in . Equation (|2~H|) therefore implies a small ratio of topological defects to pinned 
vortices, N d{ <C iVp in , for weak pinning forces compared to the elastic forces, / pin <C fi b. Sub- 
stitution of the estimate for the shear modulus quoted earlier—^ yields the field dependence 



N d{ /N pin = B cp /B for this ratio ©, where B cp = (y3/2)[7r/ln(n df aS)- 1 ] 2 (4/ p /7rJ ) 2 $ is 



the threshold magnetic field above which collective pinning holds. Last, the effect of sub- 
strate pinning by the XY model grid (jlfij) . which tends to suppress the number of unbound 
dislocations even further, is neglected here. This is valid in the regime of dilute vortex 
lattices compared to the model grid—, / < 1/36. 

Three important conclusions can be reached from the estimate for the degree of macro- 
scopic phase coherence encoded by Eqs. fj2"7j) and (J29j) above. First, observe that iV vx , A^p in 
and iVdf are all extensive thermodynamic variables that scale with the area of each layer, 
L 2 . The macroscopic phase rigidity (f2*Tj) hence attains its maximum value Jo near zero 
temperature when the total number of unbound dislocations is subthermodynamic: e.g. if 
iVdf oc L, or if A^f remains finite as L —>■ oo. This state is then a Bragg glass^^. The 
variational result ()29|) for the number of unbound dislocations obtained by Mullock and 
Evetts 2 indicates that it exists only in the absence of bulk point pins. Other types of pin- 
ning, such as surface barriers or planar defects, must therefore be present in order to impede 
flux flow by the 2D Bragg glass 16 . Second, recall that n p \ n ■ vri? 2 in ~ 1. Expression (|2*7j) 
therefore implies that weaker macroscopic phase coherence exists near zero temperature at 
dilute concentrations of unbound dislocations: p^ D \0+) > for N di <C N pin . Such a state 
is then a hexatic vortex glass 2 ^^. The variational result (J2*9~)) for the number of unbound 
dislocations indicates that this state exists at weak pinning, / pin fi b, which occurs at 
large magnetic fields B ^> B cp . Third, expression ()27j) also implies that a pinned vortex 
liquid that shows no macroscopic phase coherence is possible in the zero-temperature limit 
at sufficiently high concentrations of unbound dislocations: p^ D \0+) = if N di ~ N pin . 
This phase is then a (pinned) hexatic vortex liquid^^. The variational result ()29j) for the 




(29) 
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number of dislocations indicates that such a phase-incoherent state can occur in the regime 
of strong pinning forces, at low fields compared to the collective pinning threshold B cp . We 
remind the reader that expression (|27|) was derived by neglecting the interactions in between 
dislocations. This approximation, and Eq. (J27|) as a result, may not necessarily be valid in 
the regime of relatively dense dislocations last discussed. 

V. DISCUSSION 

The results of the previous sections are summarized by Table H] Here the low-temperature 
phases are listed by increasing order of the random pinning force, / pin . Below, we confront 
Table IJ with previous theoretical work on 2D vortex matter. 

Order Parameters. It is natural to ask what order parameters characterize the phases 
listed in Table |U Let us begin by defining the auto-correlation function 



for the Bragg-glass order parameter, where the overbar denotes a bulk average over f'. 
Drawing the analogy with the thermal degradation of phase coherence in the pristine 2D 
vortex lattice^ implies that GbtG^T) is n °t short range in the zero-temperature limit if no 
unbound dislocations are quenched in. The latter is consistent with with the properties of 
the Bragg glass listed in Table |U (see ref.— ). We conclude that the Bragg-glass phase studied 
here displays conventional phase coherence at long range. Again in analogy with the case 
of thermal disordering of the pristine vortex lattice^, the unbound dislocations that are 
quenched inside of the hexatic vortex glass, on the other hand, will result in short-range 
order in GBrdf) over a scale set by the density of such defects, n df . The absence of Bragg- 
glass order in the hexatic vortex glass, as defined above by Eq. (J3*U|) . is then consistent with 
the presence of unbound dislocations in the vortex lattice^. 

Following Fisher, Fisher and Huse 7 , we can next define the vortex glass auto-correlation 
function 



In the zero-temperature limit, this function is notably identical to unity if the groundstate 
configuration, e J< ^°, is unique. Recall now that the present treatment of the hexatic vortex 
glass has been restricted to the limit of vanishing dislocation density, in which case a unique 




dl-A/a]) 



(30) 



G VG (f) = |(expz[0(r') - 0(f + r')]>| 2 . 



(31) 
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groundstate is to be expected. The auto-correlation function GVc(r) then is not short 
range for the case of the hexatic vortex glass in the limit of weak disorder pinning. The 
zero-temperature configuration may not be unique, on the other hand, for the case of the 
(pinned) hexatic vortex liquid. The number of unbound dislocations is comparable to the 
number of pinned vortices in such case. Interactions in between dislocations may become 
important, and these may frustrate the confining action on the dislocations by the pinned 
vortices. The end result could be multiple groundstates that lead to a vortex glass auto- 
correlation function Gvg{^) that shows only short-range order. 

Consider again the (pinned) hexatic vortex liquid state that is possible at relatively dense 
concentrations of unbound dislocations, N di ~ N pin , by expression (J27)l for the superfluid 
density. Some fraction of the pairs of five-fold and seven-fold coordinated disclinations that 
make up the unbound dislocations present in this state will unbind in the strong pinning 
limit, /pi n ^> fiob. This is confirmed by direct Monte Carlo simulations^ of the Coulomb 
gas ensemble (|T7j) . Such a state then shows only short-range translational and orientational 
order. It also cannot have a net superfluid density, pf "* > 0, since it is yet more disordered 
than the hexatic vortex liquid from whence it originates. For the same reason, it can 
neither show long-range Bragg glass nor vortex glass order in the respective auto-correlation 
functions (|30J) and (|31|). This strongly pinned state is then a "conventional" vortex liquid 
(see Table |TJ). 

Vortex Glass in 2D? Fisher, Fisher and Huse argue in ref.- that the vortex-glass state 
is not possible in two dimensions above zero temperature. This statement conflicts at first 
sight with the phase-coherent vortex lattice state with a dilute concentration of quenched- 
in dislocations that we have discovered in the uniformly frustrated XY model ()16|) in the 
limit of weak random pinning. This state also shows vortex-glass order (}3Tj) ! Study of 
ref.— reveals that the authors presume the strong-pinning limit, however. The vortex glass 
is necessarily amorphous under such conditions, where it possesses a net concentration of 
unbound disclinations. It hence lies in the same topological class as the "conventional" 
vortex liquid state discussed above and listed in TablelU No conflict then truly exists between 
ref.— and the present results concerning the impossibility of observing an amorphous vortex 
glass in two dimensions. The latter is topologically distinct from the hexatic vortex glass 
discovered here. 

The hexatic vortex glass will also have some concentration of bound pairs of dislocations 
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quenched into the vortex lattice in the zero-temperature limit. Vinokur and co-workers have 
argued that the absence of infinite potential barriers along the corresponding glide planes 
will result in thermally activated plastic creep of magnetic flux due to the diffusion of such 
pairs of dislocations, and in thermally activated electrical resistance as a result^i. Hence, 
although a stack of uncoupled sheets of hexatic vortex glass shows magnetic screening in 
direct proportion to the superfluid density, pf D \ this system may not in fact be a perfect 
conductor! Flux creep immediately becomes neutralized, however, once macroscopically big 
layers (compared to the Josephson penetration depth) are coupled through the Josephson 
effect. This is due to the binding of pairs of dislocations into "quartets" that carry no 
net magnetic flux^i. The possibility just raised of three-dimensional vortex matter that is 
ohmic, but that nevertheless shows macroscopic phase coherence, therefore remains unreal- 
istic. 



VI. CONCLUSIONS 



We have demonstrated that incompressible vortex lattices that are confined to two di- 
mensions and that are void of unbound disclinations show macroscopic phase coherence near 
zero temperature in the limit of weak random pinning. The latter ensures that the total 
number of unbound dislocations quenched into the vortex lattice by the randomly pinned 
vortices is small in comparison to the total number of such pins. This in turn results in 
a net superfluid density. The hexatic vortex glass predicted here is consistent with the 
observation of isolated dislocations that are quenched into the vortex lattice of extremely 
layered high-temperature superconductors at low temperature^, and with the observation 
of superconductivity in 2D Josephson junction arrays in external magnetic field^. How 
exactly such a superconducting vortex lattice transits into a vortex liquid with increasing 
temperature remains unclear. Equations and (J2*T|) indicate that the superfluid density 
vanishes either once the quenched-in dislocations delocalize and begin to cross the length 
of the vortex lattice, or once thermally activated pairs of dislocations unbind and begin to 
cross the length of the vortex lattice^. Both mechanisms cause plastic creep of the 2D 
vortex lattice^, which destroys macroscopic phase coherence. Continuity implies that the 
melting temperature Tj 2D ) of the defective vortex lattice lies near that of the pristine vortex 
lattice^, kBT^ D "> = Jo/20, in the limit of weak random pinning. 
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APPENDIX: PINNED ELASTIC STRING 

Consider a tense elastic string of length L that lies along the x axis under periodic 
boundary conditions, and suppose that only transverse displacements, u(x), along the y 
axis are allowed. Suppose further that the string is pinned down at N P [ n sites: 

u(xi)=Vi for i = 1,2, A/p in . (A.l) 

The shape of the string with the lowest energy can then be determined by minimizing the 
elastic energy along with appropriate terms that enforce the constraints: 

1 rL / d>>i \ 2 fL 

E = -/j dx\—\ +/ dx\-[u-v]. (A.2) 

2 Jo \ax J Jo 

Here /xo denotes the shear modulus, while the field X(x) = J2i=i Xi5(x—Xi) is weighted by the 
Lagrange multipliers Aj that correspond to each of the constraints (jA.l|) . The configuration 
that minimizes the elastic energy (|A.2|) satisfies the field equation 



d 2 u 

~^ + A = (A.3) 



everywhere. Equation (jA.2|) can be easily minimized in the wave representation, 

u (x) = Y. u i eiqx and x (%) = EV i<?x > ( A - 4 ) 

i q 

expressed as a sum over allowed wavenumbers q that are multiples of ±2n/L. One then 
obtains the solution 

AW AW 

u{x) = Y.Y.G{x-x i )GT]v v (A.5) 

i=l j=l 

where G(x) = ^ q e iqx jLq 2 is the Greens functions in one dimension, and where G~j is 
the inverse of the A^n x Af p i n matrix G{x{ — x 3 ). Notice that the above solution satisfies 
the constraints (jA.l|) . It also satisfies the field equation (jA.3|) . with Lagrange multipliers 
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Aj = —/J'oY^f=iGi,j v j that weight the field X(x) at each of the pins located at X{. Direct 
substitution of the solution (|A.5|) into the elastic energy (|A.2|) then yields the result 

TV ; N in 

S = (^/2)EE«* G! y«i- ( A - 6 ) 
i=l j=l 

It reduces to the expression 

£ 1>2 = (/i /2)K - u 2 )7ki - sal (A.7) 

for the elastic energy in the special case of two pins at the thermodynamic limit, L — > oo. 
Here we used the result = G(0) — (|x|/2) for the Greens function in one dimension, 

with a limiting value G(0) —>■ oo for the constant. 
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disorder index 


phase 


p<? D) (0+)/J 


unbound dislocations? 


unbound disclinations? 


1 


Bragg glass 


unity 


no 


no 


2 


hexatic vortex glass 


fraction 


yes 


no 


3 


hexatic vortex liquid 


zero 


yes 


no 


4 


vortex liquid 


zero 


yes 


yes 



TABLE I: Listed are physical properties that characterize the low-temperature phases of vortex 
matter in two dimensions at the extreme type-II limit. The list is indexed by increasing levels of 
disorder. 
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